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Abstract 



The decay rates of cc and bb through two-photon or two-gluon annihi- 
lations are obtained by using totally relativistic decay amplitudes and a so- 
phisticated quantum-chromodynamic potential model for heavy quarkonia. 
Our results for the photonic and gluonic widths of the ^o, 3 Po) an d the 3 i"2 
states are in excellent agreement with the available experimental data. The 
procedures and mathematical techniques used by us for the treatment of the 

fermion-antifermion bound states are also applicable to other decay processes. 
12.39.Pn,13.20.Gd,13.25.-k,13.40.Hq 
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I. INTRODUCTION 



The decay rates of heavy quarkonia through two-photon or two-gluon annihilations were 
first obtained in the nonrelativistic approximation and found to be inadequate in the 
light of the experimental data ||. Improvements to the earlier results, therefore, have been 
explored by various authors by including the relativistic corrections 0-0. The quarkonium 
decay is usually treated by a suitable adaptation of the matrix element for the annihilation 
of a free quark-antiquark pair, and this treatment is either supported by appealing to the 
Bethe-Salpeter approach with instantaneous approximation |5j or simply regarded as an 
artifice @. 

In view of the ambiguity involved in the treatment of the quarkonium decays, we shall 
obtain the decay rates by following two different approaches based on different assumptions 
regarding the role of the potential energy in a bound state. In both approaches, fully 
relativistic matrix elements for the quark-antiquark annihilation will be used. Moreover, 
unlike the earlier authors, we shall use a realistic quarkonium potential model || which has 
proved highly successful in the investigation of the cc and bb spectra with spin splittings. 
Our results for the decays of the ^So, 3 Pq, and 3 P2 states of cc and bb will be compared with 
the earlier results of other authors as well as with the available experimental data. 

The general procedure and its applications to the S and P states of the fermion- 
antifermion bound states are described in Sec. II and III, while the quarkonium photonic 
and gluonic widths according to two different approaches are obtained in Sees. IV and V, 
which is followed by a discussion of our results in Sec. VI. 

II. FERMION-ANTIFERMION BOUND-STATE DECAYS 

Let us first consider the annihilation of a pair of an electron and a positron of four- 
momenta p and q into two photons of four-momenta k and k! . The second-order contribution 
of the scattering operator for this process is 
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It is to be noted that 
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It is possible to convert the matrix element from the Dirac form to Pauli form without 
making any approximation by the substitutions || 
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together with the charge-conjugation relation 

a* T ■ 

p s = a s ia 2 . 

Then, after reducing the products of the a matrices, it is found that 

F = a^ia 2 Oa r 

with 
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where we have used the center-of-mass relations 
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q = -p, g = Po, k' = -k, k' = k , (2.10) 

as well as 

k = A; k, k-e = k-e' = 0, (2.11) 
but avoided the use of the energy conservation relation 

Po = k . (2.12) 
Now, let ^ denote a positronium wave function, which can be Fourier decomposed as 

#(r) = - — / rfp^(p)e ip - x . (2.13) 
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Following the usual approach, we assume that the decay amplitude for the p component of 
the positronium wave function can be obtained from the free electron-positron annihilation 
amplitude ( |2.8|) by ignoring the relation ( |2.12| ), and treating p as a variable given by 

po = (m 2 + p 2 ) 1/2 , 0<|p|<oo, (2.14) 

while 

1 

ko = -jTYlee- ( 2 - 15 ) 

This assumption implies that the effect of the bound-state potential energy is simply to 
nullify the energy conservation relation for the free-state annihilation amplitude. 
Consequently, the amplitude for positronium decay into two photons is 
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and the resulting decay rate is given by 



dpF^(p), (2.16) 



r(«-77) = ^/^ k |l^. (2.17) 

Furthermore, the decay rates for quarkonia can be obtained by using the quarkonium 
wave functions, setting ko = ^Mqq, and making the usual multiplicative replacements in 
( [2. 17] ). F° r the decay rate T(QQ — > 77), the replacement is 



a 1 -> Ne%a\ N = 3, (2.18) 



while T(QQ — > gg) can be obtained from T(QQ — > 77) by the replacement 



e%a 2 ^ 2 -a 2 s . (2.19) 



III. S AND P STATE DECAY RATES 

We shall apply the treatment of Sec. II to obtain the decay rates for those S and P 
states of positronium which can decay into two photons, and for this purpose we shall use 
the wave functions in the matrix representation, given in Appendix A. 

A. ^Sq decay 



For the 1 Sq state ( |A10| ), the decay amplitude, given by ( |2.16| ) and ( |2.8|) , takes the form 
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With the substitution of ( |3.2| ) in ( |2.17| ), we have 
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and, upon summation over the final polarization states, 



k ■ e x e' 
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while, in view of the indistinguishability of the two photons, 

J dfi k = 2tt. 

Thus, the decay rate is given by 
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This agrees with the result in Ref. 0, where the authors have obtained the decay rate for 
the state but not for the 3 P and 3 P 2 states. 



B. decay 

For the 3 Pq state ( |A12| ), the decay amplitude ( |2.16| ) becomes 
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where only terms linear in cr in O contribute to the trace. After trace evaluation and 
simplification, we find 
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and, upon angular integrations with the help of 
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■/u po 
and Ai and A 2 are given by (|B4j). 

Substituting ( |3.10j ) in ( j2.17| ), and summing over the final polarization states, we arrive 
at the decay rate 
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C. 3 P 2 decay 



For the 3 P 2 state ( |A14| ), the decay amplitude fl2.16| ) is 

3 X 1 / 2 



where again only terms linear in cr in O contribute to the trace. 

After trace evaluation, angular integrations with the help of ( |BT 
applications of the relations 

= o, ( e , e ; - = 0, 



(3.13) 



D, and (gD, and 



(3.14) 



it is found that 



it! ( e i e 'i + e 'i e i)h + ( e ■ e')hkjh 



(3.15) 



where 



f oo 

h = I dp 



poo 

I 4 = dp 



-pM 



2p 



Po(m + p ) 



B, 



0(p), 



{A\ + A 2 ) + 9/ koP - (fr + 5B 2 + 2B V 



2p 



Po 



pl{m + p ) 



Po{m + po) 



-Bn 



(3.16) 



<f>(p), (3.17) 



7 



and Ai and Bi are given by ( p4|) and (p5|). 

Furthermore, upon averaging over the initial states with different values of M by using 
( |A8|) , and summing over the final polarization states, we obtain 

\^/'^^\h? + \h-h\% (3.18) 

which, when substituted in Q2.17D , gives the decay rate 

a 2 

T( 3 P 2 - 77) = (6|/ 3 | 2 + \h ~ h\ 2 )- (3.19) 



IV. QUARKONIUM PHOTONIC AND GLUONIC WIDTHS 

The quarkonium photonic and gluonic widths, which are obtainable by making the 
replacements ( p.!8| ) and ( |2.19| ) in the results of Sec. Ill, are given by 
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where 

^0 = \m qQ . (4.2) 

We have computed these widths by using the wave functions and parameters obtained 
from our quantum-chromodynamic potential model for heavy quarkonia ||. An essential 
feature of our model is the inclusion of the one-loop radiative corrections in the quantum- 
chromodynamic potential, which is known to be responsible for the remarkable agreement 
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between the theoretical and experimental results for the spin splittings in the cc and bb 
spectra. Another advantage of our model is that it is based on a nonsingular form of the 
quarkonium potential, and thus avoids the use of an illegitimate perturbative treatment. 

In addition to the wave functions, the parameters used for the computation of the 
widths are 

a s {cc) = 0.316, m c = 2.088 GeV, 

M(t] c ) = 2.979 GeV, M( Xc o) = 3.415 GeV, M(x c2 ) = 3.556 GeV, (4.3) 

and 

a s (bb) = 0.283, m b = 5.496 GeV, 

M(rtb) = 9.417 GeV, M( X w) = 9.860 GeV, M(x b2 ) = 9.913 GeV, (4.4) 
where we have included our theoretical value for the mass of the unobserved energy level rjf,. 



Our results for cc and bb, together with the available experimental data ]T0,|TT], are given in 
Table I. 



V. ALTERNATIVE TREATMENT OF BOUND-STATE DECAYS 

Finally, we shall explore an alternative treatment of the bound-state decays by making 
an assumption regarding the role of the potential energy which differs from that in Sec. II. 
Let us again consider the Fourier decomposition 

*(') = (ajp / dp^(p)e^ (5.1) 

of a bound-state \I/(r), and look upon it as a superposition of plane waves of different 
momenta but the same energy. Such a viewpoint is appropriate for a bound state because a 
wave packet consisting of waves of the same energy does not spread in time. It also allows 
us to treat the decay of the p component of a bound state of cc or bb as the annihilation of 
a pair of free quark and antiquark of momenta p and — p, whose energy and effective mass 
are 
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Po = \m qQ} (5.2) 



™ = (p§-P 2 ) 1/2 =(>Sq-P 2 ) • (5.3) 



and 

A 

This approach implies that the quark and antiquark in a quarkonium can be looked upon 
as free particles of constant energy and variable momentum and mass. 

The above treatment also leads to the quarkonium decay rates given by ([4.1|) in Sec. IV, 
but with an important difference. In Sec. IV, m is a constant, while p^ is a variable, given 
by ( |2.14| ). On the other hand, in the alternative treatment po is a constant, while m is a 



variable, and they are given by ( |5T2"| ) and ( |5.3| ). The computed photonic and gluonic widths 
resulting from the alternative treatment are also given in Table I. 

VI. DISCUSSION 

We have obtained the two-photon and the two-gluon relativistic decay rates of cc and 
bb by using two different approaches, which are based on apparently reasonable but very 
different assumptions. It is interesting to find that both approaches give quite similar results. 
As shown in Table I, our results for the 3 Po and 3 P2 states are in agreement with the Particle 
Data Group |T(J , while our results for the two-gluon decay of the state disagree with the 



Particle Data Group but agree with the more recent findings of the E760 collaboration ||TT 

In Table I, we have also included the nonrelativistic results obtained in Ref. [H] with the 
use of the Cornell potential. The nonrelativistic decay rates are much smaller than the ex- 
perimental values, and this disagreement has not been resolved by the authors in Refs. 
who found that the relativistic corrections amount to a reduction in the nonrelativistic decay 
rates. 

Our treatment differs from those of the earlier authors in several respects: 

1. We have used totally relativistic decay amplitudes instead of making nonrelativistic 
approximations or retaining only the leading relativistic corrections. 
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2. We have used a sophisticated quarkonium potential instead of simpler potentials such 
as the harmonic oscillator or the Cornell potential. 

3. We have used a nonperturbative treatment for the spin-dependent interaction terms in 
the quarkonium potential instead of obtaining their contributions to the energy levels 
through first-order perturbation. 



It is interesting that our values of the heavy quark masses, given by (|4.3|) and Q4.4|) , 



are somewhat higher than those generally found in the literature. This is a consequence 
of the fact that nonperturbative treatment of the spin-dependent interaction terms in the 
quarkonium potential yields constituent quark masses which are higher than those resulting 



from the commonly used perturbative treatment |T2| . The nonperturbative treatment also 
has a pronounced effect on the quarkonium wave functions. 

The remarkable agreement between our theoretical results and the experimental data 
represents a distinct success of our relativistic treatment of the bound-state decays as well 
as of our quantum-chromodynamic quarkonium model. The procedures and mathematical 
techniques used in this paper can also be applied to other decay processes. 
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APPENDIX A: FERMION-ANTIFERMION WAVE FUNCTIONS 

The fermion-antifermion wave functions in the matrix representation are given by 

*0Sb) = f^°*Kr)> ( A1 ) 
*( 3 Si) = ^itfW(r), (A2) 
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where £ 4 M is a unit vector, and is a symmetric and traceless unit tensor, such that 



E C M* C M r 
Si s, = Oij, 



M=-l 



and 



M* M 1 1 - 

M=-2 



(A3) 
(A4) 
(A5) 
(A6) 

(A7) 

(A8) 



It is also to be noted that since we treat the spin-dependent interaction terms nonperturba- 
tively, the radial wave function 0(r) has a different form for each of the above states. 

It is straightforward to show that these wave functions have the desired quantum num- 
bers by applying the operators Lj, Sf, and Jf = Lf + Sf + 2LiSi to them, and keeping in 

mind that in the matrix representation of the wavefunctions 
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In the momentum space, the corresponding wave functions are 
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APPENDIX B: ANGULAR INTEGRATIONS OF DECAY AMPLITUDES 



Angular integrations of complicated integrals appearing in the decay amplitudes can be 
performed by setting 
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TABLES 

TABLE I. Photonic and gluonic widths of cc and bb. The first two sets of theoretical results 
correspond to the relativistic treatments of Sec. IV and Sec. V. We also give the nonrelativistic 
results from Ref. Q. The experimental results for XcO an d Xc2 ar e from Ref. [10|, and those for rj c 



are from Ref. [11 1. 
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